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THERMODYNAMIC  PROPERTIES  OF  HELIUM  FOR  TEMPERATURES  BETWEEN  0° 
AND  150°  C  AND  PRESSURES  TO  300  ATMOSPHERES  GENERATED  FROM 
THE  PVT  DATA  ON  HELIUM  OBTAINED  BY  MICHELS  AND  WOUTERS 
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n  ,  .  1/ 

Dalton- 


and  Robert  E 


w  .  2/ 

Bar ieau— 


ABSTRACT 

The  isotherms  of  helium  published  by  Michels  and  Wouters  have 
been  reevaluated  for  an  equation  of  the  form 


where  the  best  values  for  the  seven  parameters  were  taken  to  be  the 
least  squares  values. 

From  our  retreatment  of  Michels  and  Wouters'  data,  we  generated 
six  ..isenthalpic  curves  along  with  compressibility  factors,  Joule- 
Thomson  coefficients,  relative  molal  heat  capacities,  relative  molal 
enthalpies,  and  molal  heat  capacities. 

1/  Research  Chemist,  Helium  Research  Center,  Bureau  of  Mines, 
Amarillo,  Texas. 

2/  Supervisory  Research  Chemist,  Project  Leader,  Thermodynamics, 
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Work  on  manuscript  completed  October  1965. 


..  .  -  •  '  '• 

'  '  *  *• ' 


1 

.  . 


-  :-pa  .  m  -  -r 


i  ■  ■  ■■  ■ 

. 


5 


INTRODUCTION 

The  Helium  Research  Center  is  presently  engaged  in  a  critical 

examination  of  all  of  the  PVT  data  on  helium  that  appear  in  the 

literature.  The  long-range  objective  is  the  development  of  a 

single  equation  of  state  for  helium  that  will  allow  all  of  the 

thermodynamic  properties  to  be  calculated  and  will  reproduce  all 

of  the  data  within  the  accuracy  with  which  the  data  are  known.  A 

preliminary  objective  is  to  obtain  the  best  values  of  second  virial 

coefficients  and  to  develop  an  equation  that  will  reproduce  the 

best  values  as  a  function  of  temperature. 

3/ 

In  a  previous  article  (,3)—  ,  we  evaluated  second  and  third 


3/  Underlined  numbers  in  parentheses  refer  to  items  in  the  list 


of  references  at  the  end  of  this  report 


virial  coefficients  of  helium  from  the  equation 

(PV)  =  A  +  B(T)P  +  C  (T)  P2 

corr.  t 


(1) 


We  then  represented  these  B's  by  a  temperature  function  similar  to 
that  employed  by  Keyes  (4) . 


B(T) 


(2) 


The  B's  evaluated  from  equation  (2)  represented  the  original  B 
values  within  0.2%  or  better.  This  was  interpreted  to  mean  that 
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equation  (2)  represented  the  second  virials  of  helium  for  the  tem¬ 
perature  range  0°  to  150°  C. 

We  have  previously  shown  (.3)  that  the  C's  can  be  represented 
by  the  equation 


C  (T) 


c . 
i 


,1/4 


+ 


+ 


3/4  m5/4 


(3) 


and  that  C's  evaluated  from  equation  (3)  are  within  the  limits  of 
about  one  standard  deviation  of  the  original  C  values.  This  was 
interpreted  to  mean  that  equation  (3)  represented  third  virials 
within  the  apparent  precision  of  the  original  data  from  which  these 
values  were  determined. 

Therefore,  since  we  had  what  we  considered  to  be  adequate  tem¬ 
perature  functions  for  representing  second  and  third  virials  of 
helium,  we  decided  to  reevaluate  all  of  the  original  data  reported 
by  Michels  and  Wouters  (6) .  An  equation  of  the  form 


was  applied  to  the  complete  set  of  data  where  the  seven  unknowns: 

(A  /T  )  ;  b  •  b  •  b  •  c.;  c  •  c  are  t.o  be  determined.  The  best 
oo  iZoiza 


values  for  the  parameters  appearing  in  equation  (4)  were  taken  to 


be  the  least  squares  values,  each  observation  being  assigned  unit, 
weight.  Our  retreatment.  of  the  original  PVT  data  reported  by 
Michels  and  Wouters  is  given  below. 
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RECALCULATION  OF  THE  ORIGINAL  HELIUM  ISOTHERMS 
REPORTED  BY  MICHELS  AND  WOUTERS 

The  experimental  data  as  published  by  Michels  and  Wouters  are 
given  in  table  1  of  this  report,  the  PV  factors  being  expressed  in 
Amagat  units  and  the  pressures  in  atmospheres  [Amagat  units  are  eval¬ 
uated  from  the  ratio  pv/p  v  ,  where  p  and  v  are  the  pressure  and 

o  o 

volume,  respectively,  at  the  high  pressure  and  at  some  given  tempera¬ 
ture  and  p  and  v  are  the  corresponding  quantities  for  the  same  mass 
o  o 

of  gas  at  0°  C  and  1  atmosphere] .  An  equation  of  the  form 


was  applied  to  the  complete  set  of  data  where  the  seven  unknowns  of 
equation  (4)  were  determined.  The  best  values  for  the  parameters 
appearing  in  equation  (4) ,  and  given  in  table  2  of  this  report,  were 
taken  to  be  the.  least  squares  values,  each  observation  being  assigned 
unit  weight.  The  deviations  given  in  table  2  are.  standard  deviations. 
These  deviations  were  determined  by  expressing  each  of  the  constants 
in  terms  of  the  original  observations  and  applying  the  law  for  the 
"Propagation  of  Errors"  (.2,  5)  . 

If  the  data  reported  by  Michels  and  Wouters  contained  a  syste¬ 
matic  error,  one  would  not.  expect  PV  to  reduce,  to  the  reference  state 
of  the  gas  [the  reference  state  is  defined  to  be:  PV  =  1  at  0°  C 
and  P  =  1  atmosphere]  .  For  example,  suppose  the  high  pressures  were, 
too  high  by  some  factor,  say  two.  Then  PV  would  reduce  to  two  at 


. 

s 


.  ' 
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TABLE  1 .  -  Experimental  PVT  data  for  helium  obtained  by  Michels 

and  Wouters 


Pressure 

atm 


9.0948 

11.2729 

14.0727 

17.6397 

21.8017 

26.7198 

31.5272 

35.1200 

36.2641 

42.6534 

52.3285 

64.2578 

78.6351 

98.010 

121.759 

152.160 

189.813 


9.9266 
12.3036 
15.3591 
19.2526 
23.7946 
29.1618 
34.4061 
38.3280 
39.5749 
46 . 5446 
57.0978 
70.1080 
85.7848 
106.923 
132.802 
165,938 
206.963 


(PV) 

exp 

Michels  and  Wouters 
Amagat.  units 


Temp.  =  0°  C 

1.00431 
1.00547 
1.00693 
1.00876 
1.01094 
1.01351 
1.01606 
1.01790 
1 . 0 1855 
1,02183 
1.02685 
1.03303 
1.04045 
1,05032 
1.06245 
1,07775 
1.09665 

Temp.  =  25°  C 

1 .09616 
1,09740 
1.09898 
1.10100 
1 . 10335 
1 .10614 
1 , 10884 
1.11088 
1.11154 
1.11505 
1.12044 
1,12708 
1.13505 
1,14584 
1 . 15881 
1 . 17533 
1,19573 


(PV)  -  (PV) 
exp 

HRC 

Amagat  units 


0.000041 
0.000057 
0.000047 
0.000006 
0.000005 
0.000002 
0.000041 
0.000007 
0.000060 
0.0000.13 
0.000008 
0.000014 
0.000023 
-0.000040 
-0.000003 
-0,0000^7 
0  000057 


0.000039 
0.000045 
0.000040 
0,000042 
0,000040 
0.000056 
0.000049 
0.000067 
0,000084 
0,000008 
-0,000019 
-0.000035 
-0.000053 
0.000023 
-0  000037 
-0,000059 
0,000079 


1/  (PV)  was  evaluated  from  equation  (4) ,  using  the  least  squares 

O  3  J.  •  #  *  -i  _  j 

values  for  the  constants  given,  in  table  2  of  this  report. 


. 
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TABLE  1.  -  Experimental  PVT  data  for  helium  obtained  by  Michels 

and  Wouters  (Con.) 


Pressure 

atm 


10.7582 
13.3340 
16.6452 
20.8640 
25.7858 
31.6005 
37.2837 
41.5319 
42.8832 
50.4317 
61.8667 
75.9532 
92 .9345 
115.820 
143.833 
179.686 
224.074 


(PV) 

exp 

Michels  and  Wouters 
Amagat  units 


Temp.  =  50°  C 

1.18799 

1.18930 

1.19100 

1.19315 

1.19568 

1.19864 

1.20158 

1.20374 

1.20446 

1.20817 

1.21402 

1.22105 

1.22965 

1.24118 

1.25507 

1.27271 

1.29459 


(PV) 


exp 

HRC 


(PV) 


ca  1 


Amagat  units 


0.000019 
0.000006 
0.000007 
-0.000006 
0.000003 
-0.000010 
0.000028 
0.000022 
0.000053 
-0.000079 
-0.000034 
-0.000132 
-0.000089 
-0.000031 
-0 . 000087 
-0.000143 
0.000066 


11.5902 

14.3650 

17.9319 

22.4766 

27.7776 

34.0412 

40.1626 

44.7368 

46.1926 

54.3225 

66.6350 

81.8034 

100.084 

124,722 

154.876 

193.458 

241.128 


Temp.  =  75°  C 

1.27986 
1.28126 
1.28307 
1.28537 
1.28804 
1.29122 
1,29436 
1 . 29663 
1.29741 
1.30138 
1.30759 
1.31510 
1 . 32425 
1.33658 
1.35143 
1.37025 
1.39312 


0.000042 

0.000030 

0,000028 

0.000021 

0.000003 

0.000012 

0.000058 

0.000019 

0,000065 

-0.000060 

-0.000036 

-0.000121 

-0.000083 

0,000033 

0.000041 

0,000043 

-0.000073 


1/  (PV)  was  evaluated  from  equation  (4),  using  the  least  squares 

value!  for  the  constants  given  in  table  2  of  this  report. 


. 
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TABLE  1 «  -  Experimental  PVT  data  for  helium  obtained  by  Michels 

and  Wouters  (Con.) 


Pressure 

atm 


12.4218 
15  .3958 
19.2184 
24.0889 
29.7694 
36.4817 
43.0393 
47.9414 
49.4995 
58.2133 
71.4029 
87.6561 
107.240 
133.625 
165.899 
207.183 
258.187 


(PV) 

exp 

Michels  and  Wouters 
Amagat  units 


Temp .  =  100°  C 

1.37170 

1.37320 

1.37512 

1.37757 

1.38040 

1.38379 

1.38707 

1.38951 

1.39029 

1.39459 

1.40115 

1.40919 

1.41893 

1.43199 

1.44761 

1.46747 

1.49168 


L(PV)exp  "  (W)cal. 
HR.C 

Amagat  units 


0.000037 

0.000038 

0.000034 

0.000034 

0.000011 

0.000035 

0.000032 

0.000022 

0.000023 

-0.000023 

-0.000024 

-0.000038 

0.000044 

0.000168 

0.000087 

0.000058 

-0.000001 


13.2532 

16.4261 

20.5043 

25.7004 

31.7595 

38.9193 

45.9135 

51.1443 

52.8042 

62.0999 

76.1779 

93.5026 

114.390 

142.512 

176.899 

220.871 

275.189 


Temp.  =  125°  C 

1.46350 

1.46510 

1.46713 

1.46973 

1.47268 

1.47625 

1.47970 

1.48234 

1.48311 

1.48770 

1.49485 

1.50318 

1.51354 

1.52723 

1.54360 

1.56442 

1.58991 


-0.000004 
0.000010 
0.000004 
0.000014 
-0.000053 
-0  000041 
-0.000061 
-0.000011 
-0.000063 
-0.000066 
0.000150 
-0.000018 
0.000149 
0.000205 
0 . 000044 
-0.000049 
-0.000043 


JL /  (PV)  was  evaluated  from  equation  (4) ,  using  the  least  squares 

C  3  x.  • 

values  for  the  constants  given  in  table  2  of  this  report. 


■ 

11 


TABLE  1 .  -  Experimental  PVT  data  for  helium  obtained  by  Michels 

and  Wouters  (Con.) 


Pressure 

atm 


14.0847 

17.4566 

21.7902 

27.3116 

33.7509 

41.3568 

48.7883 

54.3454 

56.1089 

65.9811 

80.9284 

99.3435 

121.520 

151.375 

187.909 

234.562 

292.153 


(PV) 

exp 

Michels  and  Wouters 
Amagat.  units 


Temp .  =  150°  C 


1.55532 

1.55701 

1.55914 

1.56187 

1.56502 

1.56871 

1.57235 

1.57512 

1.57593 

1.58068 

1.58807 

1.59708 

1.60787 

1.62221 

1.63967 

1.66139 

1.68792 


(PV) 


exp 


HRC 


Amagat  units 


-0.000023 
-0.000004 
-0.000020 
-0.000020 
-0.000049 
-0.000107 
-0.000121 
-0.000079 
-0  000133 
-0.000216 
-0.000119 
-0.000051 
0.000030 
0,000058 
0.000169 
-0  000005 
-0.000095 


1/ 


(PV)  1  was  evaluated  from  equation  (4) ,  using  the  least  squares 

CHI  . 

values  for  the  constants  given  in  table  2  of  this  report:. 


' 

; 
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TABLE  2 .  -  Virial  coefficients  of  helium  as  evaluated  by  the  Helium 
Research  Center  from  the  original  PVT  data  obtained  by 
Michels  and  Wouters:  0°  C  =  273,15°  K 


(A  /T  ) 
o  o 


3 


A 


o 


2 


=  (0,36590932  +  0,00000411)  x  10~2 

~3 

=  (3 , 2 .1079494  ±  0  16219855)  x  10  ' 

=  - (2 . 76069068  +  0.60075141)  x  10~ 

=  (1.64599720  +  0.55397920)  x  10”1 

=  0  99948130  +  0.00001122 

=  -(3,47366275  +  0.89251908)  x  10" 

=  (1 . 19284347  ±  0.33335061)  x  10 “4 

=  - (1 . 11341304  +  0,31003730)  x  10" 


t,°  C  Atl/ 


(-C 


10 


7,3/ 


0 

0 . 999481 

± 

0  ,000011 

0,52714 

+ 

0,0006.1 

0,8177 

+ 

0.0377 

25 

1.090959 

+ 

0.000012 

0  o  52078 

+ 

0.00042 

0,7216 

+ 

0  0226 

50 

1 . 182436 

± 

0.000013 

0.51521 

+ 

0,00041 

0,6688 

+ 

0,0195 

75 

1.273913 

+ 

0.000014 

0.51024 

+ 

0,00040 

0 , 6455 

+ 

0.0182 

100 

1.365391 

+ 

0.000015 

0,50573 

+ 

0,00038 

0 . 6426 

+ 

0.0158 

125 

1.456868 

+ 

0.000016 

0.50161 

+ 

0,00038 

0,6538 

+ 

0.0143 

150 

1,548345 

+ 

0.000017 

0,49779 

± 

0,00045 

0,6750 

± 

0.0176 

1/ 

A  =  (A  /T  )T,  where  T 
t  x  o  o' 

=  (273,1: 

5  +  t,  °  C)  , 

2/  Evaluated  from  equation  (2) ,  using  the  least  squares  values  for 
b^,  ,  and  b^  given  above, 


3/  Evaluated  from  equation  (3) ,  using  the  least  squares  values  for 
c^y  c^,  and  c^  given  above. 


f 
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P  =  1  atm.  and  0°  C.  Therefore,  whether  or  not  our  function,  equa¬ 
tion  (4) ,  reduces  to  the  reference  state  of  the  gas  can  be  used  to 
test  for  the  presence  of  a  systematic  error. 

We  carry  out  this  test,  for  a  systematic  error  in  the  following 
way.  We  evaluate  the  unknowns  of  equation  (4)  by  least  squares 
methods.  We  then  evaluate  our  function,  equation  (4),  at  T  =  273.15°  K 
and  P  =  1  atm.  and  compare  this  value  to  one.  If  the PV  product,  is 
different  from  one,  a  statistical  test,  is  carried  out  to  see  if  this 
difference  is  significant.  If  this  difference  is  of  statistical 
significance,  then  this  value,  represents  the  systematic  error  which 
is  removed  by  multiplying  all  of  the  published  PV!s  by  the.  appropriate 
factor  to  remove  it. 

Suppose  we  evaluate  the  right-hand  side,  of  equation  (4)  at.  0°  C 
and  P  =  1  atm.  From  the  data  given  in  table.  2,  we  have 


An  =  V  =  0.99948130  +  0,00001122,  Amagat:  units 

0°  \T  /  o 
o 


V 


"  (tTm 


+ 


+ 


2-3  = 


o 


T  3/4  '  T  5/4' 
o  O 


(0,52714  ±  0,00061)  x  10 


•3 


co° 


C1  °2 

+ 


''3—\  = 


T  1/,‘4  T  i/4  T  5/4 
o  o  o 


-  (0,8.177  +  0.0377)  x  .10 


where  T  =  273,15°  K,  Therefore, 
o 


PV 


Ao°  +  Bo°  +  co° 


o°  c 


)P=1  atm. 


1,00000836 


/ 
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Since  PV  is  not  exactly  equal  to  one,  this  indicates  that  there 
might  be  a  systematic  error  in  the  data.  Now  in  order  to  decide  if 
this  difference  is  statistically  significant,  we  need  to  have  an 
estimate  of  the  standard  error  of  PV  at  0°  C  and  P  =  1  atm.  We  eval¬ 
uated  the  standard  error  of  our  function,  equation  (4) ,  at  0°  C  and 
P  =  1  atm.  by  expressing  A^0 ,  Bq0 ,  and  C^0  in  terms  of  the  original 
data  and  applying  the  law  for  the  "Propagation  of  Errors"  (2,  5), 
When  we  did  this,  we  got 


S  ! 

(PV) 

jj0°  c 

n 

P=  1  a  t  in . 


+  0  00001095 


or 


I  0°  C 
P=  1  atm. 


1.000008  +  0.000011,  Amagat  units 


We  see  that  the  difference  between  the.  reference  state  of  the.  gas 
and  the  PV  product  evaluated  from  equation  (4)  is  less  than  a  stand¬ 
ard  deviation.  This  is  taken  as  statistical  evidence  that  there  is 
no  systematic  error  in  the  experimental  data  reported  by  Michels  and 
Wouter s . 

We  have,  previously  shown  (3)  that  the  temperatures  as  given  by 
the  primary  temperature  scale  have  smaller  uncertainties  than  those 
evaluated  from  the  PV  measurements.  We.  concluded,  therefore,  that. 


/ 


. 

r 
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the  primary  scale  was  more  accurate  and  that  the  temperature  of  each 
isotherm  is  as  given  by  the  primary  scale.  This  means  that  we  should 
set  A  =  RT  for  each  isotherm,  where  T  is  evaluated  from  the  primary 
temperature  scale. 

We  reevaluated  the  PVT  data  given  in  reference  6  in  the  following 
way:  from  equation  (4),  for  T  =  T^  =  273.15°  K  and  P  =  1  atm.,  we 
have 


or  , 


o  o  o  o  o  o 


Substituting  equation  (6)  into  equation  (4) ,  we  get 


o  o 


+  C2(t3/4  '  T  7/4)  +  C3(t5/4  '  T  9/4) 

-  o  o 


Equation  (7)  was  then  applied  to  the  complete  set  of  data  to 
give  the  residuals  given  in  column  3  of  table  3.  The  parameters 
appearing  in  equation  (7)  were  evaluated  by  least  squares  solution, 


.  .  .  .  '  .  IS":  .  •  •  "■  6 

i  ..  !j‘  "'C  '■  "  ■  V..  \  1,  •  : 

■  .  .  j  i  '  i i 

K  v  f..  ;  : 

■  j  •  "  •  ■.  .  ■ . . .  £  J  ■■■  ''  ;i  .• '  •  ■  su  > .. .  i 


■C  ‘‘ 


- 


. 


■  .  •  -  :■  .  -  •  pa 

n  u  l  >g.  ■  ,  p  '  Vi.-  • 
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TABLE  3 .  -  Experimental  compressibility  data  for  helium  as 

reevaluated  using  equations  (6)  and  (7): 

PV  =  1  at  P  =  latm.  and  0°  C 


Pressure 

atm 


9.0948 
11.2729 
14.0727 
17.6397 
21.8017 
26.7198 
31.5272 
35.1200 
36.2641 
42,6534 
52.3285 
644  2578 
78.6351 
98.010 
121.759 
152.160 
189.813 


9.9266 
12.3036 
15.3591 
19.2526 
23.7946 
29.1618 
34.4061 
38.3280 
39.5749 
46 . 5446 
57.0978 
70. 1080 
85.7848 
106.923 
132.802 
165.938 
206.963 


(PV) 

exp 

Michels  and  Wouters 
Amagat  units 


Temp.  =  0°  C 

1.00431 

1,00547 

1.00693 

1.00876 

1.01094 

1.01351 

1.01606 

1.01790 

1.01855 

1.02183 

1.02685 

1.03303 

1.04045 

1.05032 

1.06245 

1.07775 

1.09665 


I  (PV)  -  (PV)  . 
L  exp  cal 

HRC 

Amagat  units 


0.000048 

0.000063 

0.000052 

0.000011 

0.000010 

0.000006 

0.000044 

0.000009 

0.000062 

0.000015 

0.000008 

0.000013 

0.000021 

-0.000043 

-0.000006 

-0.000049 

0.000060 


Temp.  =  25°  C 

1.09616 
1.09740 
1.09898 
1. 10100 
1. 10335 
1.10614 
1.10884 
1.11088 
1.11154 
1.11505 
1.12044 
1.12708 
1. 13505 
1.14584 
1.15881 
1.17533 
1.19573 


0.000046 

0.000052 

0.000046 

0.000047 

0.000045 

0.000059 

0.000051 

0.000069 

0.000086 

0.000009 

-0.000019 

-0.000037 

-0.000056 

0.000019 

-0.000041 

-0.000061 

0.000083 


j  £'  '  • 

■» .  »  »  ■  ,  . 


f 

- 

■  ■ 

. 
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TABLE  3 .  -  Experimental  compressibility  data  for  helium  as 

reevaluated  using  equations  (6)  and  (7): 

PV  £  1  at  P  =  1  atm,  and  0°  C  (Con.) 


Pressure 

atm 


10.7582 
13 . 3340 
16.6452 
20.8640 
25.7858 
31.6005 
37.2837 
41.5319 
42.8832 
50.4317 
61 . 8667 
75.9532 
92.9345 
115.820 
143.833 
179.686 
224.074 


11.5902 

14.3650 

17.9319 

22.4766 

27.7776 

34.0412 

40.1626 

44.7368 

46.1926 

54.3225 

66.6350 

81.8034 

100.084 

124.722 

154.876 

193.458 

241.128 


Michels  and  Wouters 
Amagat  units 

Temp .  =  50°  C 

1.18799 

1.18930 

1.19100 

1.19315 

1.19568 

1.19864 

1.20158 

1.20374 

1.20446 

1.20817 

1.21402 

1.22105 

1.22965 

1.24118 

1.25507 

1.27271 

1.29459 

Temp.  =  75°  C 

1.27986 
1.28126 
1.28307 
1.28537 
1.28804 
1.29122 
1 . 29436 
1.29663 
1.29741 
1.30138 
1.30759 
1.31510 
1.32425 
1.33658 
1.35143 
1.37025 
1.39312 


(PV)  -  (PV)  . 

L  'exp  'cal. 

HRC 

Amagat  units 


0.000027 

0.000013 

0.000013 

0.000000 

0.000008 

-0.000006 

0,000031 

0.000024 

0.000055 

-0.000078 

-0.000035 

-0.000134 

-0.000093 

-0.000035 

-0.000091 

-0.000145 

0.000071 


0.000050 
0.000038 
0 . 000035 
0.000027 
0.000008 
0.000016 
0.000061 
0.000022 
0.000067 
-0.000059 
-0.000037 
-0.000123 
-0.000087 
0,000028 
0.000037 
0,000041 
-0.000068 


: 

*  •  •  *  • 
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TABLE  3.  -  Experimental  compressibility  data  for  helium  as 

reevaluated  using  equations  (6)  and  (7): 

PV  =  1  at  P  =  1  atm,  and  0°  C  (Con.) 


Pressure 

atm 


12.4218 

15.3958 

19.2184 

24.0889 

29.7694 

36.4817 

43.0393 

47.9414 

49.4995 

58.2133 

71.4029 

87.6561 

107.240 

133.625 

165.899 

207.183 

258.187 


13.2532 

16.4261 

20.5043 

25.7004 

31.7595 

38.9193 

45.9135 

51.1443 

52.8042 

62.0999 

76.1779 

93.5026 

114.390 

142.512 

176.899 

220.871 

275.189 


(PV) 

exp 

Michels  and  Wouters 
Amagat  units 


Temp.  =  100°  C 

1.37170 

1.37320 

1.37512 

1.37757 

1.38040 

1.38379 

1.38707 

1.38951 

1.39029 

1.39459 

1.40115 

1.40919 

1.41893 

1.43199 

1.44761 

1.46747 

1.49168 

Temp .  -  125°  C 

1.46350 

1.46510 

1.46713 

1.46973 

1.47268 

1.47625 

1.47970 

1.48234 

1.48311 

1.48770 

1.49485 

1.50318 

1.51354 

1.52723 

1.54360 

1.56442 

1.58991 


<W>e*p  -  (PV) 
HRC 

Amagat  units 


0.000046 

0.000047 

0.000041 

0.000041 

0.000016 

0.000039 

0.000035 

0.000024 

0.000025 

-0.000022 

-0.000025 

-0.000040 

0.000039 

0.000163 

0.000082 

0.000056 

0.000004 


0.000005 

0.000019 

0.000012 

0.000021 

-0.000047 

-0.000037 

-0.000057 

-0.000009 

-0.000060 

-0.000065 

0.000149 

-0.000021 

0.000144 

0.000200 

0.000039 

-0.000051 

-0.000038 


cal . 


• 

; 

' 

1  ... 

, 

• 

• 

S':  ■  .  ; 

■ 
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TABLE  3.  - 


Pressure 

atm 


14.0847 

17.4566 

21,7902 

27.3116 

33.7509 

41.3568 

48.7883 

54.3454 

56.1089 

65.9811 

80.9284 

99.3435 

121.520 

151.375 

187.909 

234.562 

292.153 


Experimental  compressibility  data  for  helium  as 

reevaluated  using  equations  (6)  and  (7)  : 

PV  ^  1  at.  P  =  1  atm,  and  0°  C  (Con.) 


(PV)  „ 

(PV)  -  (PV) 

exp 

L  exp 

Michels  and  Wouters 

HRC 

Amagat  units 

Amagat  unit 

Temp.  =  150°  C 

1,55532 

-0.000013 

1.55701 

0.000005 

1.55914 

-0.000012 

1.56187 

-0.000013 

1.56502 

-0.000043 

1.56871 

-0.000102 

1.57235 

-0.000117 

1.57512 

-0.000076 

1.57593 

-0.000131 

1.58068 

-0.000215 

1.58807 

-0.000120 

1.59708 

-0.000054 

1.60787 

0.000025 

1.62221 

0.000053 

1.63967 

0.000163 

1.66139 

-0.000007 

1.68792 

-0.000090 

■ 
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each  observation  being  given  unit  weight.  These  constants  are  given 
in  table  4  of  this  report.  The  deviations  given  in  table  4  are  stand¬ 
ard  deviations.  These  deviations  were  evaluated  by  expressing  each  of 
the  constants  in  terms  of  the  original  observations  and  applying  the 
law  for  the  "Propagation  of  Errors"  (2,  5) . 

From  our  reevaluation  of  Michels  and  Wouters'  data,  we  evaluated 
thermodynamic  properties  of  helium  for  temperatures  between  0°  and 
150°  C  and  pressures  to  300  atmospheres.  Our  calculations  and  the 
quantities  we  evaluated  are  given  in  the  next  section. 

THERMODYNAMIC  PROPERTIES  OF  HELIUM  AS  EVALUATED  FROM  THE  RESULTS 

OF  MICHELS  AND  WOUTERS'  WORK 

We  generated  six  isenthaipic  curves  from  the  results  of  our  recal¬ 
culation  of  Michels  and  Wouters'  PVT  data  on  helium,  corresponding  to 
the  isenthaipic  zero  pressure- temperature  intercepts  reported  by  Barieau 
and  Briggs  (^L)  in  their  reevaluation  of  Roebuck  and  Osterberg's  Joule- 
Thorns  on  data  for  helium.  We  generated,  these  isenthaipic  curves  from  Michels  and 
Wouters'  data  in  the  following  way.  The  helium  isotherms  reported  in 
reference  6  were  represented  by  the  equation 


where  A^0  is  defined  by  equation  (6)  and  the  b's  and  c's  are  the  least 
squares  values  previously  determined  from  equation  (7) . 


: 
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TABLE  4.  -  Vlrial  coefficients  for  helium  as  calculated  from 

equations  (6)  and  (7) 


b  =  (3.20923530  +  0.16188694)  x  10 

b2  =  -(2.74977122  +  0.59947520)  x  10 

b3  =  (1.63456005  +  0.55275691)  x  lo" 

C;L  =  -(3.45274355  +  0.89045538)  x  10 

c2  =  (1.18441022+0.33255397)  x  10“ 

c3  =  -(1.10625748  +  0.30932473)  x  10 


=  1 


1  atm. 


t 


T  3/4  T  5/4 
o  o 

(B  x  103)+ 


T  1/4  T  3/4  T  5/4 
o  o  o 

-(C  x  101)- 


0 

0.999473 

25 

1.090949 

50 

1.182426 

75 

1.273902 

100 

1.365379 

125 

1.456855 

150 

1.548332 

1/  B  (T) 

1 

ml/4  L 

0.52735  +  0.00054 
0.52101  +  0.00031 
0.51544  +  0.00028 
0.51047  +  0.00026 
0.50597  +  0.00022 
0.50184  +  0.00022 
0.49802  +  0.00033 

A.  +  !ai 

t1/2  Tj 


0.8273  +  0.0354 
0.7310  +  0.0189 
0.6778  +  0.0154 
0.6541  +  0.0142 
0.6507  +  0.0117 
0.6614  +  0.0104 
0.6820  +  0.0150 


2/  C(T) 


The  variances  and  covariances  of  the  parameters  appearing  in 
equation  (7)  are  given  in  the  appendix  of  this  report. 


Note : 


.  >  ,  ■■  :  ' 


- 


. 
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Dividing  both  sides  of  equation  (8)  by  TP, 


V 


=  /V\  /_n_  + 

\T  P )  +  0/4 


r'-.  +  ^9/4)  +  (^74  +  p/4  +  p/4)1,  (9) 


Differentiating  equation  (9)  with  respect  to  (1/T) ,  keeping  the 
pressure  constant,  we  get 


d(V/T) 

Is  ( 1/T) 


a] 

9J 


4T 


1/4 


7b„  9b. 


5bX  + 


,1/2  T 


9b3\  / 

T  )  +  (5C1 


7c,  9c,, 

+  ~U2  +—' |P 


(10) 


But 


5  (V/T) 


d(l/T) 


— *P 


■  ( 


BH\ 

SP/r 


(11) 


Substituting  equation  (11)  into  equation  (10)  ,  we  get 


4T 


1/4 


(5bi 


7b  9b 

i  — ± - 1-  — - 

ml/2  T 


)  +  (5ci 


7c„  9c 


+  ,1/2  +  T  ^ 


ft 


(12) 


Integrating  equation  (12)  from  P  -  0  to  P  -  P  at  constant  temper/ 
ture ,  we  get 


( 


HP  “  HP=0 


4T 


1/4 


(5bi 


7b  9b 

.  —4—  +  — - 

_  1  /  2  T 


)p  +  (be 


7c  9c  n  2 
2  ,  „3\  P_ 

1/2  T  /  2 


(13) 


We  employed  equation  (13)  and  the  expression 


(AH)p  =  C°  AT  =  — 


2 


R  AT 


(14) 


/ 
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to  generate  isenthalpic  curves  for  helium  from  Michels  and  Wouters’ 
data  corresponding  to  the  isenthalpic  zero  pressure- temperature 
intercepts  reported  in  reference  1.  The  isenthalpic  zero  pressure- 
temperature  intercepts  reported  by  Barieau  and  Briggs  (_1)  are  given 
in  table  5  of  this  report;  the  isenthalpic  curves  we  generated  from 
the  data  of  Michels  and  Wouters,  corresponding  to  the  zero  pressure- 
temperature  intercepts  of  table  5,  are  given  in  table  6  of  this  re¬ 
port. 

Compressibility  factors  for  helium  for  temperatures  between  0° 
and  150°  C  and  pressures  to  300  atmospheres  are  given  in  table  7  of 
this  report.  These  Z' s  were  evaluated  by  dividing  both  sides  of 
equation  (8)  by  the  factor  (Aq0/T^)T,  When  we  do  this,  we  see  that 
Z  is  expressible  by  the  equation 


im 


(V/To 


i  + 


i 

<W 


+ 


(15) 


The  deviations  listed  in  table  7  are  standard  deviations .  These 
deviations  were  determined  by  expressing  Z  in  terms  of  the  original 
observations  and  applying  the  law  for  the  "Propagation  of  Errors 
(2,  5) . 

Table  8  of  this  report,  gives  pR  products  as  a  function  of  tem¬ 
perature  and  pressure,  where  p  is  the  density  and  R  is  the  gas  con 
stant.  These  pR  products  were  determined  from  the  equation 


pR 


JL 

ZT 


(16) 


. 


. 
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TABLE  5.  -  Isenthalpic  zero  pressure- temperature  intercepts  reported 

bv  Briggs  and  Barieau  in  their  recalculation  of 

Roebuck  and  Osterberg's  lx  data  for  helium 

for  temperatures  between  0°  and  150  C 


Zero  Pressure-Temperature  Intercept 

„  o  r  °  K 

Curve  c 

438.564 
386.249 
361.364 
336.939 
309.920 
285.291 


150 

165.414 

100 

113.099 

75 

88.214 

50 

63.789 

25 

36.770 

n 

12.141 

. 


- 


T 

°K 

438, 

438 

437 

436 

435 

434 

433 

432 

431 

430 

429 

428 

427 

426 

425 

424 

423 

422 


Isenthalpic  curves  generated  from  the  PVT  data  on  helium  obtained  by 


Michels 

and  Wouters 

p 

T 

P 

T 

P 

atm 

°  K 

atm 

°K 

atm 

0.0000 

386. 

249 

0.0000 

361.364 

0.0000 

9.2506 

386 

4.0150 

361 

5.8132 

25.6527 

385 

20.1462 

360 

21.7957 

42.0560 

384 

36.2888 

359 

37.7973 

58.4609 

383 

52.4438 

358 

53.8193 

74.8678 

382 

68.6122 

357 

69.8631 

91.2774 

381 

84.7951 

356 

85.9301 

107.6903 

380 

100.9934 

355 

102.0217 

124.1069 

379 

117.2084 

354 

118.1395 

140.5279 

378 

133.4411 

353 

134.2850 

156.9539 

377 

149.6927 

352 

150.4599 

173.3854 

376 

165.9644 

351 

166.6659 

189.8233 

375 

182.2575 

350 

182.9046 

206.2681 

374 

198.5733 

349 

199.1780 

222.7204 

373 

214.9131 

348 

215.4878 

239.1811 

372 

231.2783 

347 

231.8361 

255.6508 

371 

247.6704 

346 

248.2249 

272.1303 

370 

264.0907 

345 

264.6563 

288.6203 

369 

280.5409 

344 

281.1324 

305.1217 

368 

297.0224 

343 

297.6557 

• 

... 

V 

> 

• 
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TABLE  6. 


Isenthalpic  curves  generated 

Michels  and 


T 

P 

T 

°K 

atm 

°K 

336.939 

0.0000 

309.920 

336 

14.8420 

309 

335 

30.6751 

308 

334 

46.5379 

307 

333 

62.4320 

306 

332 

78.3594 

305 

331 

94.3221 

304 

330 

110.3221 

303 

329 

126.3616 

302 

328 

142.4427 

301 

327 

158.5679 

300 

326 

174.7395 

299 

325 

190.9600 

298 

324 

207.2321 

297 

323 

223.5585 

296 

322 

239.9421 

295 

321 

256.3858 

294 

320 

272.8927 

293 

319 

289.4662 

292 

318 

306.1096 

291 

from  the  PVT  data  on  helium  obtained  by 


Wouters  (Con,) 


P 

atm 

T 

°  K 

P 

atm 

0.0000 

285.291 

0.0000 

14.3429 

285 

4.4646 

29.9722 

284 

19.8419 

45.6449 

283 

35.2764 

61.3635 

282 

50.7715 

77.1307 

281 

66.3309 

92.9493 

280 

81.9584 

108.8224 

279 

97.6582 

124.7530 

278 

113.4345 

140.7444 

277 

129.2919 

156.8001 

276 

145.2353 

172.9234 

275 

161.2697 

189.1184 

274 

177.4005 

205.3889 

273 

193.6335 

221.7390 

272 

209.9747 

238.1732 

271 

226.4307 

254.6961 

270 

243.0083 

271.3126 

269 

259.7150 

288.0278 

268 

276.5585 

304.8472 

267 

293.5472 

ho 


- 

■ 

• 

■ 

TABLE  7 .  -  Compressibility  factors  of  helium;  Z  =  P/pRT 


T,  ° 

K 

P  ,atm. 

0 

273  . 

15 

1.0000000 

275 

1*0000000 

300 

1*0000000 

325 

1*0000000 

350 

1*0000000 

375 

1.0000000 

400 

1.0000000 

425 

1.0000000 

P ,  a  tm . 
1 


1.0005275 

+ 

0.0000005 

1.0005235 

+ 

0.0000005 

1.0004742 

+ 

0.0000003 

1.0004331 

+ 

0.0000002 

1.0003983 

+ 

0.0000002 

1.0003685 

+ 

0.0000002 

1.0003426 

+ 

0.0000002 

1.0003200 

+ 

0.0000002 

P  ,  a  tm . 
5 


1 .  002636 1 

+ 

0.0000026 

1.0026159 

+ 

0.0000025 

1.0023695 

+ 

0.0000013 

1.0021641 

+ 

0.0000011 

1.0019903 

+ 

0.0000010 

1.0018414 

+ 

0.0000008 

1.0017122 

+ 

0.0000008 

1.0015993 

+ 

0.0000011 

T,  ° 

K 

TABLE  7.  -  ' 

Compress ibility 

factors  of  helium:  Z  -  P/pRT 

(Con. ) 

P  ,atm. 

50 

P  ,  a  tm . 

10 

P 

,atm. 

20 

-J 

273  • 

15 

1,0052680 

+ 

0.0000051 

1,0105194  + 

0,0000095 

1.0261744 

+ 

0,0000190 

275 

1,0052277 

+ 

0.0000048 

1.0104390  + 

0.0000089 

1.0259756 

+ 

0.0000178 

300 

1,0047357 

+ 

0 , 0000026 

1.0094581  + 

0.0000048 

1.0235462 

+ 

0.0000098 

325 

1,0043254 

+ 

0,0000022 

1.0086395  + 

0,0000042 

1.0215137 

+ 

0.0000088 

350 

1,0039781 

+ 

0,0000019 

1,0079461  + 

0,0000036 

1.0197887 

+ 

0.0000076 

375 

1.0036804 

+ 

0.0000015 

1.0073512  + 

0,0000028 

1.0183069 

+ 

0.0000060 

400 

1.0034222 

+ 

0.0000015 

1.0068354  + 

0.0000028 

1,0170206 

+ 

0,0000060 

425 

1,0031964 

+ 

0.0000021 

1.0063839  + 

0.0000040 

1.0158939 

+ 

0.0000086 

TABLE  7.  - 

C  ompr  essibility 

factors  of  helium:  Z  -  P/pRT 

(Con. ) 

T,  °  K 

P  ,  a  tm . 

75 

P  ,  a  tm . 

QQ 

P ,  a  tm . 

150 

273  o 1 5 

1 *0391064 

+ 

0.0000229 

1.0519349  ± 

0.0000242 

1.0772815  ± 

0.0000285 

275 

1.0388110 

+ 

0.0000215 

1.0515446  ± 

0.0000227 

1.0767070  ± 

0.0000267 

300 

1 o  035 1953 

+ 

0.0000121 

1.0467618  + 

0.0000129 

1.0696468  ± 

0.0000142 

325 

1.0321641 

+ 

0.0000111 

1 .0427435  ± 

0.0000122 

1.0636893  ± 

0.0000130 

350 

1.0295874 

+ 

0.0000096 

1.0393224  + 

1 

0.0000107 

1.0586010  ± 

0.0000114 

375 

1.0273713 

+ 

0.0000076 

1.0363765  + 

0.0000086 

1.0542089  ± 

0.0000092 

400 

1.0254460 

+ 

0.0000079 

1.0338148  + 

0.0000090 

1.0503826  ± 

0.0000098 

425 

1.0237584 

+ 

0.0000113 

1.0315679  + 

0.0000129 

1.0470221  + 

0.0000141 

TABLE  7.  -  Compressibility 


P  ,  a  tm . 

T,  °  K  200 


273*15 

1  *  1022142 

+  0*0000554 

275 

1 « 1014628 

+  0*0000520 

300 

1  *  09220 1 1 

+  0*0000247 

325 

1  ©  0843513 

+  0*0000183 

350 

1  *  0776246 

±  0*0000158 

375 

1  *  071 8041 

+ 0*0000122 

400 

1 *  0667241 

+ 0*0000109 

425 

1*0622564 

+  0*0000154 

factors  of  helium:  Z  =  P./pRT  (Con.) 


P  ,  a  tm .  P ,  a  tm . 

250  300 


1*1267331 

+ 

0*0001061 

1*1508381 

+ 

0*0001765 

1.1258120 

4“ 

0*0000998 

1*1497545 

± 

0.0001662 

1*1144249 

+ 

0*0000471 

1.1363180 

+ 

0*0000793 

1* 1047295 

+ 

0*0000332 

1*1248237 

+ 

0*0000563 

1*0963932 

± 

0*0000282 

1.1149067 

+ 

0.0000476 

1*0891620 

± 

0*0000211 

1*1062827 

,;+ 

0*0000356 

1*0828392 

± 

0*0000163 

1*0987280 

+ 

0*0000271 

1*0772708 

+ 

0*0000224 

1*0920654 

+ 

0.0000370 

TABLE  8.  -  Values  of 


T,  °  K 

P  ,  a  tm . 

0 

P  ,  a  tm . 

1 

273  . 15 

0 

0  *  00365906 

275 

0 

0  <>00363446 

300 

0 

0.00333175 

325 

0 

0 .00307559 

350 

0 

0.00285601 

375 

0 

0  *00266568 

400 

0 

0.00249914 

425 

0 

0.00235219 

pR ,  atm/°  K 


P ,  a  tm . 

5 

P  ,  a  tm . 

10 

0.01825683 

0.03641807 

0o 01 81 3438 

0.03617453 

0o  01662727 

0.03317622 

0.01535139 

0.03063671 

0.01425734 

0.02845822 

0.01330883 

0.02656888 

0.01247863 

0.02491474 

0.01174592 

0.02345444 

Lo 


TABLE  8.  -  Values  of  pR ,  atm/°  K  (Con.) 


T,  0  K 

P  ,atm. 

20 

P ,  a  tm . 

50 

P ,  a  tm . 

75 

P ,  a  tm . 

100 

273  *  1 5 

0  o  07245763 

0. 17838061 

0e 26424091 

0.34802460 

275 

0.07197591 

0o 17721491 

0* 26253792 

0.34581163 

300 

0o06604203 

Go  16283258 

0o2415Q032 

0.31844240 

325 

0.06101135 

0o 15060606 

0.  22357805 

0.29507959 

350 

0.05669237 

0o 14008504 

0.20812775 

0.27490439 

375 

0.05294413 

0c 13093630 

0.19467158 

0.25730675 

400 

0o04966055 

0.12290803 

0.18284727 

0.24182281 

425 

0.04676031 

0ol 1580645 

0.17237523 

0,22809368 

•  - 

4 


TABLE  8.  -  Values  of  pR ,  at.m/°  K  (Con.) 


T  0  K 

P ,  a  tm . 

?  50 

P  ,atm. 

200 

P  ,  a  tm  . 

250 

P  ,  a  tm . 

300 

273.15 

0  . 50975425 

0.  66429774 

0. 81230242 

0.95434591 

275 

0.50659513 

0.66027896 

0.80749799 

0.94881913 

300 

0.46744404 

Go  61038819 

0. 74776985 

0.88003530 

325 

0.43390343 

0.56751404 

0.69630693 

0.82064145 

350 

0  o  40484698 

0o 53026681 

0  c  65148682 

0.76880230 

375 

0c 37943144 

0o 49760339 

0.61209138 

0.72314248 

400 

0o3570 1276 

Oo 46872476 

0.57718634 

0.68260755 

425 

0o 33709049 

0.  44300816 

0.54604216 

0  c  64637368 

u> 

GJ 
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where  P  is  the  pressure  in  atmospheres,  T  is  the  temperature  in  °  K, 
and  Z  is  the  compressibility  factor  of.  helium  corresponding  to  a 
given  P  and  T. 

If  we  divide  both  sides  of  equation  (13)  by  R,  where  R  =  (An0/T  )  , 

u  o 

then  we  get  an  expression  from  which  we  can  evaluate  relative  enthal¬ 
pies,  L/R.  That  is, 


4(AQO/To)T1/4 


7b  9b  , 

5bi  +7m  +  — ,p 


7c2  9c3\  P2 

,5ci  +  371  +  — )  1 

L_  T 


4 


(17) 


Relative  enthalpies  are  given  in  table  9  of  this  report  for  tempera¬ 
tures  between  0°  and  150°  C  and  pressures  to  300  atmospheres  along 
with  the  standard  deviation  of  each  of  these  values. 

If  we  differentiate  equation  (17)  with  regard  to  the  temperature, 
keeping  the  pressure,  constant,  we.  get  an  expression  of  the  form 


d<HP  -  HP=o)/dT 

~R 


16(A0°/To)T5/4 


21b 


T 


2 
.1/2 


+ 


21c  ^ 

JJ2  * 


45c3 


(18) 


T 


'V 


TABLE  9. 

-  Values  of  relative 

enthalpies , 

P 

P=0 

R 

L  O  y 

R  ’  K 

T,  ° 

K 

P  ,  a  tm . 
0 

P  ,  a  tm . 

1 

P, 

,  a  tm . 

5 

P  ,atm. 

10 

273* 

15 

0 

0*1 64 

+ 

0.002 

0*822 

+ 

0.008 

1.642 

+ 

0*016 

275 

0 

0.164 

+ 

0*002 

0*820 

+ 

0*008 

1*639 

+ 

0.015 

300 

0 

0*161 

0*001 

0*807 

+ 

0.004 

1.613 

+ 

0*008 

325 

0 

0*159 

+ 

0.000 

0*796 

+ 

0.002 

1.591 

± 

0.004 

350 

0 

0*157 

+ 

0*000 

0*786 

+ 

0*002 

1*572 

± 

0*  004 

375 

0 

0.156 

*4- 

0.001 

0*778 

+ 

0.003 

1.556 

± 

0.006 

400 

0 

0*154 

+ 

0*001 

0*771 

± 

0.004 

1*542 

± 

0*007 

425 

G 

0.153 

+ 

0*001 

0.765 

+ 

0.004 

1.530 

+ 

0.009 

TABLE 

9. 

-  Values 

of  relative 

enthalp: 

T,  °  K 

P  ,atm 
20 

■ 

P  ,  a  tm . 

50 

273.15 

3®  277 

+ 

0.031 

8.148  ± 

0,070 

275 

3.273 

+ 

0.030 

8.139  ± 

0®  067 

300 

3.221 

+ 

0.016 

8.021  ± 

0.037 

325 

3.179 

“1- 

0,008 

7.923  + 

0.019 

350 

3.142 

± 

0.008 

7.839  + 

0.018 

375 

3.111 

+ 

GoGll 

7,766  + 

0.025 

400 

3.084 

± 

0.014 

7.702  + 

0®033 

425 

3.060 

+ 

0.017 

7.643  + 

0.039 

VHP=0  _L_ 

R  R  ’ 


P ,  a  tm . 
75 


12.166 

+ 

0.097 

12.153 

+ 

0.093 

11.992 

+ 

0.051 

11.856 

± 

0.026 

11.739 

+ 

0.025 

11.635 

+ 

G®036 

11.542 

+ 

0.046 

11.459 

+ 

0.054 

(Gon. ) 

P ,atm. 

100 

16,147  ±  0,118 
16,132  ±  0,114 
15,937  ±  0,062 
15,770  ±  0,032 
15,624  ±  0,031 
15,494  ±  0,044 
15,376  ±  0.057 
15.269  +  0.067 


U> 


TABLE  9.  -  Values  of  relative  enthalpies 


P  ,  a  tm . 


T,  ° 

K 

150 

200 

273* 

15 

23®  997 

+ 

0c  147 

31*698 

+ 

0  o  16 1 

275 

23*979 

+ 

0*  141 

31*681 

+ 

0*1  54 

300 

23o  746 

+ 

0*076 

31*449 

4- 

0*082 

325 

2  3*  540 

+ 

0®  038 

31*233 

+ 

0  ®  041 

350 

23*353 

+ 

0*039 

31*028 

+ 

0  *  044 

375 

23*183 

+ 

QoG57 

30.833 

+ 

0*065 

400 

23®  025 

+ 

0o073 

30*646 

+ 

0*083 

425 

22  o  877 

+ 

0  o  0  8  6 

30*467 

+ 

0*096 

H  H 
P  P=0 

R 


L  o 
R  ’ 


K  (Con.) 


P ,atm. 
250 


P  ,  a  tm . 
300 


39*250 
39c 237 
39*047 
38  *  849 
38*  647 
38 « 444 
38®  241 
38®039 


+  0o 168 

+  Oo 16 1 

+  0*085 
+  0o042 
±  0*046 
+  O0O08 
+  0*087 
+  0*101 


46  * 653 
46  *  647 
46*  53Y 
46  *388 
46®  211 
46*016 
45c 809 
45*593 


±  0*184 
±  0.177 
+  0®  095 
+  0  * 048 
+  0.048 
+  0o 07  0 
+ 0.090 
+  0* 105 


U> 
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from  which  the  relative  molai  heat  capacities  listed  in  table  10 
were  calculated.  The.  deviations  given  in  table  10  are  standard 
deviations . 

Suppose  we  rearrange  equation  (18),  substituting  Cp(He)  =  5/2  R. 
When  we  do  this,  we  get  the  following  expression: 


P 


R 


=  2.5 


1 


16(A0O/To)T 


5/4 


3* 

.  21b  45b 

(5b  +  .  ,  ’  -f  ~~ 

\  1  1/2  T 


-)p  +  ( 


21c 2  45c3\P2 

\5cl  4172  +  1 


7 


from  which  the  molai  heat  capacity  values  given  in  table  1.1  of  this 
report,  were  calculated.  The  deviations  given  in  table  11  are  stand¬ 
ard  deviations, 

Jouie-Thomson  coefficients  for  helium  at  zero  pressure  for 
temperatures  between  0°  and  150°  C  were  calculated  as  follows.  4p-Q 
is  defined  by  the  equation 


P=0 


IS). 


(BH/BP) 


P=0 


P=0 


(BH/BT) 


PJ 


P=0 


(20) 


But 


(BH/BT) 


Pj 


P=0 


00 

=  C.p  and  C  for  helium  =  2,5  R  Therefore, 


(BH/BP)r 


T. 


P=0 


(BH/BP) 


P=0 


.0 


2.5  R 


(21) 


We  have  previously  shown  that  (BH/BP)  is  expressible  by  equation 
(12).  Substituting  equation  (12)  into  equation  (21),  jj,  _  is 


(19) 


' 


. 


)  br 


TABLE  10. 

-  Values 

of 

relative 

molal  heat  capacities, 

Cp  Cp=0  J 

R  R 

T  ° 

K 

P ,  a  tm . 

0 

P  ,  a  tm « 

1 

P ,  a  tm . 

5 

P 

,  a  tm . 

10 

273  , 

15 

0 

-0,0001 

+ 

0,0000 

- 0  «  0006  + 

0,0002 

-0,0012 

±  0,0003 

275 

0 

-0,0001 

+ 

0,0000 

-0,0006  + 

0,0002 

-0,0012 

±  0,0003 

300 

0 

-0,0001 

± 

0,0000 

-0,0005  + 

0,0001 

-0,0010 

±0,0002 

325 

0 

-0,0001 

+ 

0,0000 

- 0 ®  0004  + 

0,0001 

-0,0008 

±  0,0002 

350 

0 

-QoGGOl 

+ 

0,0000 

-0,0003  ± 

0,0001 

-0,0007 

+ 0,0001 

375 

0 

-OoOOOl 

+ 

0,0000 

-0,0003  + 

0,0000 

-0,0006 

± 0,0001 

400 

0 

-OoGOGl 

-i- 

0,0000 

-0,0003  + 

0,0000 

-0,0005 

±  0,0001 

425 

0 

-0,0000 

+ 

0,0000 

-0,0002  + 

0,0000 

-0,0005 

+  0  ®  0000 

TABLE  10.  -  Values  of  relative  molal  heat  capacities. 


J 

R 


(Con. ) 


VCP=0 

R 


P ,  a  tm 

20 


P ,  a  tm . 

50 


P  ,  a  tm . 
75 


P ,  a  tm . 

100 


273*15 

-0.0023 

”j~ 

0*  0007 

-0*0052 

■4- 

0*0015 

-  0*0071 

Hb 

0*0021 

-0.0086 

+ 

0.0026 

275 

-0*0023 

+ 

0.0007 

-0*0052 

+ 

0*0015 

0 

r*~ 

0 

0 

6 

O 

1 

± 

0*0021 

-0*0085 

± 

0*0026 

300 

-  0.0019 

+ 

0o  0005 

-0*0043 

+ 

0*0011 

-  0.0059 

+ 

0.0015 

-0*0072 

+ 

0*0019 

325 

-  0*0016 

+ 

0*0003 

-0.0036 

+ 

0.0008 

-  0*0050 

± 

0*0011 

-0*0062 

0.0014 

350 

-  0.0013 

+ 

0*0003 

-0*0031 

+ 

0*0006 

-  0*0044 

0  *0008 

-0.0055 

+ 

0*0010 

375 

-  0*0012 

+ 

0.0002 

-0*0027 

+ 

0*0004 

-  0*0039 

+ 

0*0006 

-0*0049 

4- 

0*0007 

400 

-  0*0010 

+ 

0*0001 

-0*0024 

4. 

0*0003 

-  0*0035 

+ 

0*0004 

-0*0045 

+ 

0*0005 

425 

-  0*0009 

0.0001 

-0.0022 

+ 

0*0002 

-  0.0032 

+ 

0*0003 

-0*0041 

+ 

0*0004 

-P' 

o 


TABLE  10.  - 

Values  of 

relative  molal  heat  capacities, 

C 

p"cp=o 

R 

“  (Con.) 

P, 

atm. 

p 

,a  tm . 

P, 

atm. 

P 

,  a  tm . 

T  ^ 

K 

150 

200 

250 

300 

273  o 

15 

-0.0101 

+ 

0o0033 

-0  0  0096 

4 

0.0037 

-  0*  0073 

4 

0.0039 

-0.0031 

± 

0.0041 

275 

-0.0100 

+ 

0.0033 

-0o  0096 

4 

0.  0036 

-  0*0074 

4 

0.0038 

-0.0033 

+ 

0.0040 

300 

-0  . 0087 

4 

0o0024 

1 

0 

0 

O 

O 

00 

s0 

4- 

0oQ026 

-  0.0078 

4 

0.0027 

-0.0053 

4 

0.0029 

325 

-0.0078 

4 

0  0  00 17 

-0.0084 

4 

0o0019 

-  0.0080 

+ 

0.0020 

-0.0066 

+ 

0.0021 

350 

-  0  o  0071 

+ 

0o0013 

-Q.0080 

4 

0.0014 

-  0.0081 

“f- 

0.0015 

-0.0075 

4 

0.0016 

375 

-0.0066 

+ 

0o0009 

-0. 0076 

+ 

OoOOlO 

-0.0081 

4 

0.0011 

-0.0081 

4 

O.GGll 

400 

-O0OO6I 

4 

0e 0007 

-0o0073 

4 

0.0007 

-  0.0081 

+ 

0.0008 

-0.0085 

+ 

0.0008 

425 

-0o0057 

+ 

0o0005 

-0.0070 

4 

0o0005 

-0.0080 

4 

0.0005 

-0.0087 

4 

0.0006 

T,  ° 

K 

P  ,atm. 
0 

273* 

15 

2  o  5000 

275 

2  o  5000 

300 

2o 5000 

325 

2  o  5000 

350 

2 o  5000 

375 

2  o  5000 

400 

2  ® 5000 

425 

2  o  5000 

11.  -  Values  of  molal 


P  ,atm 

l 


o  4999 

+ 

OoOOOO 

®  4999 

-f 

OoOOOO 

o4999 

+ 

OoOOOO 

o  4999 

4- 

OoOOOO 

e  4999 

4 

OoOOOO 

«  4999 

4- 

OoOOOO 

o  4999 

OoOOOO 

o  5000 

+ 

OoOOOO 

TABLE 

2 

2 

2 

2 

2 

2 

2 

2 


heat  capacities,  ~ 
~  K 


P ,atm. 
5 


2c 4994 

+ 

0*0002 

2  ® 4994 

+ 

0o0Q02 

2  » 4995 

+ 

OoOOOl 

2  o  4996 

+ 

OoOOOl 

2  o  4997 

+ 

OoOOOl 

2<,4997 

4- 

OoOOOO 

2o4997 

+ 

OoOOOO 

2  o  4998 

4- 

OoOOOO 

I 

*  ,atm. 

10 

2o4988 

1  4- 

0o0003 

2*  4988 

f  + 

0*0003 

2  *  4990 

0o0002 

2  0  4992 

+ 

0o0002 

2*4993 

+ 

OoOOOl 

2o  4994 

OoOOOl 

2  0  4995 

± 

OoOOOl 

2  0  4995 

+ 

OoOOOO 

TABLE  ii„  -  Values  of  moiai  heat  capacities, 


(Con. ) 


P  ,  a  tm .  P ,  a  tm . 

20  50 


P ,  a  tm . 


P  ,atm. 
U1CL 


273,15 

2  o  4977 

+ 

0,0007 

275 

2,4977 

± 

0  o  0007 

300 

2  , 4981 

4~ 

0o0005 

325 

2  o  4984 

+ 

0,0003 

350 

2,4987 

+ 

0  o  0003 

375 

2,4988 

+ 

0o0002 

400 

2 , 4990 

4 

OoOOOl 

425 

2 ,4991 

4- 

OoOQOl 

2,4948 

+ 

0o0015 

2  o4929 

2  o4948 

4 

0c0015 

2,4930 

2,4957 

± 

OoOOll 

2 ,494 1 

2,4964 

4 

Ocooos 

2  *  4950 

2  0  4969 

4 

0 0  0006 

2  0  4956 

2  o4973 

4 

0o0004 

2,4961 

2o  4976 

4 

0o0003 

2,4965 

2 ,4978 

4 

0o0002 

2,4968 

±  0,0021 

2  o4914 

4 

0,0026 

±  0o0021 

2,4915 

4 

0,0026 

±  0,0015 

2,4928 

4 

0,0019 

±  OoOOll 

2,4938 

4 

0,0014 

4  0o00G8 

2,4945 

4 

0,0010 

4  0  0  0  0  0  6 

2,4951 

4 

0,0007 

4  0o0004 

2,4955 

4- 

0,0005 

+  0,0003 

2,4959 

4 

0,0004 

TABLE  11.  -  Values  of  molal  heat 


T,  ° 

K 

P 

,  a  tm . 

150 

P, 

atm. 

200 

273. 

15 

2.4899 

± 0.0033 

2.4904 

+ 

0.0037 

275 

2.4900 

± 0.0033 

2.4904 

+ 

0.0036 

300 

2.4913 

± 0.0024 

2.4911 

+ 

0.0026 

325 

2.4922 

+  0.0017 

2.4916 

+ 

0.0019 

350 

2.4929 

± 0.0013 

2.4920 

+ 

0.0014 

375 

2.4934 

+  0.0009 

2.4924 

+ 

0.0010 

400 

2.4939 

+  0.0007 

2.4927 

+ 

0.0007 

425 

2.4943 

4~  0.0005 

2.4930 

+ 

0.0005 

capacities , 


(Con. ) 


P  ,  a  tm . 
250 


P ,a tm. 
300 


2.4927 

+ 

0.0039 

2.4926 

+ 

0.0038 

2.4922 

+ 

0.0027 

2.4920 

+ 

0.0020 

2.4919 

+ 

0.0015 

2.4919 

+ 

0.0011 

2.4919 

+ 

0.0008 

2.4920 

+ 

0.0005 

2.4969 

+ 

0.0041 

2.4967 

+ 

0.0040 

2.4947 

+ 

0.0029 

2.4934 

+ 

0.0021 

2.4925 

+ 

0.0016 

2.4919 

+ 

0.0011 

2.4915 

+ 

0.0008 

2.4913 

+ 

0.0006 

4> 

•J^ 
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expressible  as 


4P=0 


4T1/4(2.5R) 


(5bi 


7b2  9b 

+  ~m+  — 


(22) 


where  R  =  (Aqo/Tq).  We  employed  equation  (22)  in  evaluating  the 
(j,p_Q  values  for  helium  given  in  column  4  of  table  12.  The  standard 
deviations  of  the  |j,p  values  calculated  from  equation  (22)  are 
also  given  in  this  table.  These  standard  deviations  were  evaluated 
by  expressing  M-p_Q  in  terms  of  the  original  data  and  applying  the 


law  for  the  "Propagation  of  Errors"  ( 2 ,  5)  . 


f 


. 
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TABLE  12 .  -  Joule-Thomson  coefficients  calculated  from  the  original 

PVT  data  on  helium  obtained  by  Michels  and  Wouters 


Curve 

!  Zero  Pressure- 

-Tempera 

ture  Intercept 

[( 

-p,  X 

lo2)P=o 

°  C 

0  K 

°  C 

°  K/atm. 

273.150 

0.000 

6 

.578 

+ 

0.065 

275.000 

1.850 

6 

.569 

+ 

0.063 

0 

285.291 

1/ 

12. 141 

6 

.521 

+ 

0.050 

300.000 

26.850 

6 

.460 

+ 

0.035 

25 

309.920 

1/ 

36.770 

6 

.422 

+ 

0.027 

325.000 

51.850 

6 

.369 

+ 

0.018 

50 

336.939 

1/ 

63.789 

6 

.332 

+ 

0.015 

350.000 

76.850 

6 

.294 

+ 

0.016 

75 

361.364 

1/ 

88.214 

6 

.263 

+ 

0.019 

375.000 

101.850 

6 

.229 

+ 

0.023 

100 

386.249 

1/ 

113.099 

6 

.202 

+ 

0.026 

400.000 

126.850 

6 

.  172 

+ 

0.030 

425.000 

151.850 

6 

.122 

+ 

0.035 

150 

438.564 

1/ 

165.414 

6 

.097 

+ 

0.037 

1/ 

Values  reported  in 

reference  1  for  the  average 

z< 

ero  ] 

pressure 

temperature  intercept. 
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APPENDIX 

The  variances  and  covariances  of  the  parameters  appearing  in 
equation  (7)  are 


S4 

= 

2.62073816  x 

10” 

Sb2 

= 

3.59370510  x 

10~ 

S>3 

= 

3.05540196  x 

10“ 

s2 

C1 

= 

7.92910782  x 

io‘ 

s2 

C2 

— 

1.10592146  x 

10" 

s2 

C3 

= 

9.56817913  x 

10" 

s2 

blb2 

= 

-9.69954948  x 

10 

sb  b 

= 

8.92894480  x 

10“ 

blb3 

Sb  c 

= 

-1.33199188  x 

10 

bl  1 

Sb  c 

= 

4.95015561  x 

10" 

bl°2 

Sb  c 

= 

-4.57577779  x 

10 

bl  3 

s2 

b2b3 

= 

-3.31177972  x 

10 

Sb  c 

— 

4.95085186  x 

10" 

2  1 

s2 

= 

-1.84191718  x 

10 

b2°2 

-7 


-10 


■8 


-4 


-7 


r 

49 


"v. 

=  1.70451669  x  10"6 

*Vi 

=  -4.57756903  x  10"8 

b3°2 

=  1.70495217  x  10"6 

s2 

b3C3 

=  -1.57961978  x  10"5 

s2 

C1C2 

=  -2.95963973  x  10"U 

s2 

C1C3 

=  2.74817895  x  10"10 

s2 

C2°3 

=  -1.02806132  x  10"8 

The  above  quantities  were  determined  by  expressing  each  constant  in 
terms  of  the  original  observations  and  applying  the  law  for  the 
“Propagation  of  Errors.1'  The  standard  deviation  in  any  other  cal¬ 
culated  thermodynamic  quantity  evaluated  from  equation  (8)  can  be 


determined  from  these  variances  and  covariances. 
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